Are direct photons a clean signal of a thermalized quark gluon plasma? 
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Direct photon production from a quark gluon plasma (QGP) in thermal equilibrium is studied 
directly in real time. In contrast to the usual S-matrix calculations, the real time approach is 
valid for a QGP that formed and reached LTE a short time after a collision and of finite lifetime 
(~ 10 — 20 fm/c as expected at RHIC or LHC). We point out that during such finite QGP lifetime 
the spectrum of emitted photons carries information on the initial state. There is an inherent 
ambiguity in separating the virtual from the observable photons during the transient evolution 
of the QGP. We propose a real time formulation to extract the photon yield which includes the 
CO ' initial stage of formation of the QGP parametrized by an effective time scale of formation 

This formulation coincides with the S-matrix approach in the infinite lifetime limit. It allows to 
separate the virtual cloud as well as the observable photons emitted during the pre-equilibrium stage 
from the yield during the QGP lifetime. We find that the lowest order contribution 0(ot ern ) which 
does not contribute to the S-matrix approach, is of the same order of or larger than the S-matrix 
contribution during the lifetime of the QGP for a typical formation time ~ 1 fm/c. The yield for 
momenta > 3 Gev/c features a power law fall-off ~ T 3 T 2 /k 5 and is larger than that obtained with 
ON , the S-matrix for momenta > 4 Gev/c. We provide a comprehensive numerical comparison between 

CN . the real time and S-matrix yields and study the dynamics of the build-up of the photon cloud and 

the different contributions to the radiative energy loss. The reliability of the current estimates on 
photon emission as well as theoretical uncertainties on the details of the initial state are discussed. 
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Amongst the different potential experimental signatures of the formation and evolution of a quark gluon plasma 
(QGP) that is conjectured to be formed in ultrarelativistic heavy ion collisions, hard electromagnetic probes, namely 
direct photons and dileptons are considered to be very promising pi 0- Photons and dilepton pairs only interact 
1} ■ electromagnetically and their mean free paths are much larger than the size of the QGP, thus these electromagnetic 
probes leave the hot and dense region after formation without further scattering, carrying with them clean information 
of the early stages of the collision. Therefore a substantial effort has been devoted to obtaining a theoretical assessment 
of the spectra of direct photons and dileptons emitted from a thermalized QGPpJ-Q. Preliminary assessments 
concluded that direct photon emission from a thermalized QGP can be larger than that from the hadronized phase0,0|, 
sparking an intense effort to obtain reliable estimates of the direct photon spectrum [H, 0, II] • For recent reviews on 
theoretical and phenomenological aspects of electromagnetic probes, see [lOt IllL Il2j . 

The first observation of direct photon production in ultrarelativistic heavy ion collisions has been reported by 
the WA98 collaboration in 208 Pb + 208 Pb collisions at sfs = 158 Gev at the Super Proton Synchrotron (SPS) at 
CERN 13]. The results display a clear excess of direct photons above the expected background from hadronic decays 
in the range of transverse momentum pt > 1.5 Gev/c in the most central collisions. These results provide an 
experimental confirmation of the feasibility of direct photons as reliable probes in ultrarelativistic heavy ion collisions 
and may pave the way for understanding the formation and evolution of a QGP. 

A variety of fits of theoretical results to the experimental data had been reported^!]], however, the results seem 
inconclusive, models with or without QGP emission seem to fit the data in a manner compatible with models based 
solely on hadronic 'cocktails' (for a detailed review see [ill]). 

The current ultrarelativistic heavy ion program at the Relativistic Heavy Ion Collider (RHIC-BNL) and the proposed 
heavy ion program Alice at the forthcoming Large Hadron Collider (LHC-CERN) have as a main goal to continue 
the experimental pursuit of the long-sought QGP with beam energies of y/s ~ 200 AGev for Au + Au at RHIC and 
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up to y/s ~ 5500AGev for Pb + Pb at CERN. This active experimental program with the possibility of statistical 
analysis on event-by-event basis justifies the theoretical assessment of experimental probes at a deeper level. For a 
recent summary of measurements at RHIC seefl-ij. 

The S-matrix approach to calculating the photon yield from a QGP in local thermal equilibrium treats the plasma 
as stationary and with an infinite lifetime, while it is clear that QGP is a transient, non-equilibrium state|l5l Il6|. 
Current theoretical understanding suggests that a QGP may be formed ~ 1 fm/c after a nucleus-nucleus collision 
and thcrmalizcs via parton-parton scattering. The subsequent evolution is assumed to be described by hydrodynamic 
expansion until the temperature cools down to the hadronization scale ^160 Mev. At RHIC the initial temperature 
of the plasma is expected to be of the order of 300 Mev, and assuming Bjorken's longitudinal expansion with a cooling 
law T(t) = Ti(ti/t)z it is expected that the lifetime of the QGP is of order < 10 fm/c for a hadronization temperature 
of about 160 Mev. At the LHC the initial temperature is expected to reach ~ 450 Mev and the lifetime of the QGP 
would be expected to be of order ~ 20 — 30 fm/c. The transverse size of the QGP formed in the most central collisions 
is of the order of the radius of the nucleus which for Pb + Pb is about 7 fm thus the typical space-time dimension of 
QGP in local thermal equilibrium is about 10 fm. 

Despite the fact that the quark gluon plasma, if formed, will occupy a finite and rather small volume in space time, 
the S-matrix approach to obtain the photon and dilepton yields treats the plasma as a medium in thermal equilibrium 
and of infinite extent in space-time|l5l Il6|. The production rate obtained from this approach is then input into a 
spacetime evolution combined with a hydrodynamic expansion of the plasma [ToL I lll lT ^ . A recent analysis of photon 
production along these linesQ to fit the data from WA98 13j suggests that the large px region is dominated by the 
first few fm/c of (hydrodynamic) evolution and is very sensitive to the early stages of the evolution. 

The issue of a finite space-time extension of the QGP and the hadronic phase has received attention with respect 
to the emission of photons and dileptons. The influence of a finite spatial size of the plasma has been addressed for 
the emission of thermal photons |l7l Il8| and more recently for thermal dileptons |19|| from a hadronic gas, where the 
breaking of detailed energy-momentum conservation by finite size effects was studied. 

Preliminary studies of the finite lifetime effects on the photoproduction yield were reported in rcf . 20] . The results 
of those studies pointed out the importance of non-equilibrium real time processes whose contribution is subleading in 
the infinite lifetime limit, but that are of the same order or larger than the S-matrix contribution during the lifetime 
of a QGP expected at RHIC and LHC. Two main consequences of the study in refs.[2(| are: 

i) during a finite lifetime the spectrum of direct photons is sensitive to the initial conditions that lead to a thermalized 
QGP with the large px region of the spectrum more sensitive to the initial stages, and 

ii) to lowest order a em the spectrum resulting from the non-equilibrium processes flattens for momenta px > 
2 Gev/c. The sensitivity of the large px part of the spectrum to initial conditions has also been pointed out in ref. 9], 
and perhaps coincidentally, the WA98 data^l displays a flattening of the spectrum for px > 1.5 Gev/c. 

Goals of this article: The goals of this article are to continue the study of direct photon production from a QGP 
in local thermodynamic equilibrium with a finite lifetime, directly in real time. We focus on the following aspects: 

• Assessing the contribution to the direct photon spectrum from the lowest order processes that are subleading in 
the infinite time limit. These processes are qq — > 7 and q — > (77 and correspond to the one loop contribution to 
the photon polarization, namely of order a em . The contributions of these processes vanish in the infinite time 
limit and do not contribute to the rate obtained from the S-matrix approach, but do contribute to the yield 
during a finite lifetime. Focusing on the lowest order contributions we identify the dynamical aspects of photon 
production in real time in the simplest possible case. This study highlights that there are contributions to all 
orders in a s that are being missed by the S-matrix calculation. 

• A detailed analysis of the dynamics of the build-up of the virtual photon cloud and to provide a systematic 
effective description of the initial stage between the collision and thermalization that allows a clear separation 
of the virtual photons. We discuss the inherent difficulties associated with an unambiguous identification of the 
virtual photon cloud during a finite time interval. 

• A systematic description of direct photon production during a finite time interval including the initial preparation 
of the state. 

• An analytic and numerical comparison of the real time yield obtained in lowest order, namely of 0(a em ) and 
the S-matrix yield, of order a s a em In ^-0,0 to assess the potential experimental significance of the processes 
that are missed by the S-matrix calculation but that contribute to the direct photon yield from a QGP with a 
finite lifetime. We provide a comprehensive numerical study of the direct photon yield to lowest order 0(a em ) 
with an analysis of the spectrum. 

• A study of the dependence of the spectrum on the initial conditions prior to the onset of local thermal equilibrium. 
This study reveals important aspects of the initial conditions prior to thermalization that influence the spectrum. 
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• A study of the radiative energy loss, in particular the contributions associated with the interaction energy as 
well as the cooling of the plasma by photon emission. 

• A simple energy-time uncertainty argument would suggest that for momenta larger than the inverse lifetime 
of the QGP, the effects of a finite lifetime should be subleading. Our study clearly shows this expectation not 
to bear out. In fact we show that contributions from the region uj ^ k in the imaginary part of the photon 
polarization are very important during the finite lifetime and of the same order (or larger ) than the usual result 
valid solely for lu — k even for photons with large transverse momentum. 

This article is organized as follows: in section[H]we revisit the S-matrix approach to highlight its caveats. In section 
IIIII wc present the real time formulation to photon production beginning with a full gauge invariant treatment of the 
electromagnetic interaction of quarks. In section Hvl we provide a simple and transparent derivation of the expression 
for the photon production yield in real time to lowest order in a em and finite QGP lifetime. This formulation 
reproduces the results obtained in ref . |20j by a more general kinetic description and is explicitly shown to coincide 
with the S-matrix formulation in the infinite QGP lifetime limit. In this section we address the issue of initial 
conditions and in particular the subtle but important aspects associated with the formation of the photon cloud. In 
this section we present a detailed analysis of the radiative energy loss and the different contributions, providing an 
analytic and numerical study of the total energy radiated during the lifetime of the QGP. In section[V]we address the 
issue of the electromagnetic dressing of the initial state (density matrix) by providing an initial density matrix that 
includes the photon cloud parametrized by a formation time scale of the QGP after the parton stage following an 
ultrarelativistic heavy ion collision, conjectured to be ~ 1 fm/c. This parametrization interpolates smoothly between 
the adiabatic preparation of asymptotic states and the uncorrelated initial state assumed in the S-matrix calculation. 
The consideration of such initial states (density matrix) allows us to address the issue of the formation time and 
includes in a phenomenological manner the photon cloud of the pre-equilibrium stage. 

Our conclusions are presented in section I*VT1 



In order to highlight the shortcomings of the S-matrix approach to calculate photon emission, and to establish 
contact with the real-time approach to photon production introduced in Sec. IV, we now summarize some important 
aspects of the scenario of QGP formation and evolution and the S-matrix approach to the calculation of photon 
emission. 

As mentioned in the introduction QGP is conjectured to be formed in ultrarelativistic heavy ion collisions from the 
deconfinement of strongly interacting quarks and gluons in the incoming nuclei. The details of the dynamics of the 
collision are not completely understood, nor, in particular, the electromagnetic aspects of the parton distribution 
functions. It is conjectured that immediately after the collision the partons are almost free and parton-parton 
scattering leads to a state of (local) thermal equilibrium on a time scale ~ 1 fm/c after the collision. The photons 
emitted during the pre-equilibrium stage are assumed to leave the medium |24j|. The thermalized QGP undergoes 
adiabatic hydrodynamic expansion during a lifetime of ~ 10 — 20 fm/c after which the plasma hadronizes. The QGP 
in local thermal equilibrium under the strong interactions is not in equilibrium under the electromagnetic interactions 
resulting in photons emitted directly from the thermalized plasma. 

The S-matrix approach to the calculation of photon emission begins by writing the Hamiltonian in the form 



where Hqqb is the full QCD Hamiltonian, H~ t is the free photon Hamiltonian, and Hi nt is the interaction Hamiltonian 
between quarks and photons with the quark electromagnetic current, A M the photon field, and e the electromagnetic 
coupling constant. 

Consider that at some initial time U the state \i) is an eigenstate of Hq with no photons. The transition amplitude 
at time tf to a final state |/, "f\{p)) = \ f) <8> \"fx(p)), again an eigenstate of Hq but with one photon of momentum p 
and polarization A, is up to an overall phase given by 



II. S-MATRIX APPROACH AND ITS CAVEATS 
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(II.1) 



S(t f ,t i ) = (f, 7x (f)\U(t f ,t i )\i) , 



(II.2) 
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where U(tf,ti) is the time evolution operator in the interaction representation 



-i / H intJ (t) ,dt 



U(t f) ti) = Texp 

L Jtt 

= 1 - ie J t 1 dt J d * x Ji(x,t) Af,j(x,t) + 0(e 2 ) , (II.3) 

where the subscript / stands for the interaction representation in terms of Hq. In the above expression we have 
approximated U(tf,ti) to first order in e, since we are interested in obtaining the probability of photon production to 
lowest order in the electromagnetic interaction. The usual S-matrix element for the transition is obtained from the 
transition amplitude S(tf,U) above in the limits U — > — oo and tf —* +oo 

S fi = S(+oo, -oo) = f d 3 x dt e lP " x » e x u (f\J"(x)\i) + 0(e 2 ) , (II.4) 



'2E 

where E = \p\ and = (E,p) are the energy and four-momentum of the photon, respectively, and e x is its polarization 
four-vector. Since the states \i) and |/) are eigenstates of the full QCD Hamiltonian Hqcb, the above S-matrix element 
is obtained to lowest order in the electromagnetic interaction, but to all orders in the strong interaction. We note that 
the S'-matrix element in effect is the amplitude for the transition between asymptotic states |z;in) — > |/, 7x(p); out), 
i.e., Sfi = (/, 7a(p); out|i; in), where |/, 7>(p); out) = |/;out) (g> \j\(p);out). Here, \f\(p);out) is the asymptotic out 
state with one photon of polarization A and momentum p, and \i; in) (|/; out)) is the asymptotic in (out) state of the 
quarks and gluons. 

The rate of photon production per unit volume from a QGP in thermal equilibrium at temperature T is obtained 
by squaring the S'-matrix element, summing over the final states, and averaging over the initial states with the 
thermal weight e~ l3Ei /Z(j3), where f3 — 1/T, Ei is the eigenvalue of Hq corresponding to the eigenstate \i), and 
Z((3) = J2i e ~ l3Ei 1S the partition function. Using the resolution of identity 1 = J2f l/K/l) the sum of final states 
leads to the electromagnetic current correlation function. Upon using the translational invariance of this correlation 
function, the two space-time integrals lead to energy-momentum conservation multiplied by the space-time volume 
£1 = V(tf — ti) from the product of Dirac delta functions. The term tf — ti — * +oo is the usual interpretation of 
2nS(0) in the square of the energy conserving delta functions. 

These steps lead to the following result for the photon production rate in the S-matrix approach 0, El 

dN 1 1 d 3 p ^ -0E ilQ |2 2 ixv w < , m d 3 p 

d^ = nz(p)Wffr x lM 9 w ^ [p) m^f' ( } 

where W^ v (K) is the Fourier transform of the thermal expectation value of the current correlation function defined 
by 

W<(K) = J d 4 xe lKx (J^Q)J v {x))p . (11.6) 

In the expression above (• • • ) | g denotes the thermal expectation value. To lowest order in e 2 but to all orders in the 
strong interactions, W< V (K) is related to the retarded photon self-energy II* (K) by 3 

, lmU^(uj = k,k) , 



Thus, one obtains the (Lorentz boost) invariant photon production rate 

(II. 



dN _ g» v lmR^ v {u) = k,k) 



d 3 pd 4 x (2tt) 3 e k / T -l 

All the calculations of the photon production yield from a thermalized QGP in equilibrium begin by obtaining 
lmll^ u (u> — fc, A:) to calculate the rate. The most recent result up to leading logarithmic order in the strong coupling 
has been obtained in ref.[8|. 

We have reproduced the steps leading to Eq. I|II.8H , which is the expression for the photon production rate used in 
all S-matrix calculations in the literature, to highlight several important steps in its derivation in order to compare 
and contrast to the real-time analysis discussed below. The main features of the above result that will be compared 
to the real time computation are the following: 
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• The initial states \i) are averaged with the thermal probability distribution at the initial time ti for quarks and 
gluons. In the usual calculation this initial time ti —* — oo, as emphasized above and the initial state describes 
the photon vacuum and a thermal ensemble of quarks and gluons. Thus the quarks and gluons are assumed 
to have thcrmalized in the infinite past. Furthermore, this treatment also assumes that the quarks and gluons 
are asymptotic states in the infinite past, thus neglecting the fact that these arc confined in the colliding nuclei 
before the collision. 

• The transition amplitude is obtained via the time evolution operator U(tr,ti) evolved up to a time tf and the 
transition amplitude is obtained by projecting onto a state |/) at time tf, which in the calculation is taken 
tf — > +00. The sum over the final states leads to the electromagnetic current correlation function averaged over 
the initial states with the Boltzmann probability distribution, i.e., the thermal expectation value of the current 
correlation function. 

• Taking tf — ► +00 and ti — ► — 00 and squaring the transition amplitude leads to energy conservation and an 
overall factor tf — tj. The rate (transition probability per unit time per unit volume V) is finally obtained 
by dividing by (tf — U)V . The important point here is that taking the limit of tf — ti — ► +00 results in two 
important aspects: energy conservation and an overall factor of the time interval tf — ti. The resulting rate is 
independent of the time interval and only depends on the photon energy (and obviously the temperature). 

Main assumptions in the S-matrix approach:. In order to compare our methods and results with those 
obtained within the usual S-matrix framework described above, it is important to highlight the main assumptions 
that are implicit in all previous calculations of photon production from a thcrmalized QGP and that are explicitly 
displayed by the derivation above. 

• The initial state at U (which in the usual calculation is taken to —00) is taken to be a thermal equilibrium 
ensemble of quarks and gluons but the vacuum state for the physical transverse photons. 

• The usual calculation of the rate to lowest order in the electromagnetic coupling, entails that there are no 
electromagnetic corrections to the intermediate states, namely, there is no photon dressing of the states that 
enter in the thermal density matrix. 

• Taking ti — > — 00, tf — > +00 manifestly assumes that quarks and gluons are asymptotic states in the infinite 
past and in the infinite future. Obviously this is inconsistent with the fact that before the collision quarks and 
gluons should be described in terms of their parton distribution functions in the nuclei. Furthermore assuming 
quarks and gluons to be asymptotic states as tf — > +00 manifestly ignores the hadronization phase transition 
to a confined phase at a finite time of order 10 — 20 fm/c. In ref.|9( confinement in the initial state had been 
encoded in a 'confinement factor', namely a phenomenological parameter included to account for the effects of 
confinement. 

• Assuming the QGP to have equilibrated at ti — > —00 and taking tf — > +00 makes explicit that the plasma is 
assumed to be described as a stationary state in thermal equilibrium at all times. 

• The buildup of population of photons is neglected along with the electromagnetic dressing of quarks, these as- 
sumptions are generally invoked to justify a calculation of the yield or rate to lowest order in the electromagnetic 
coupling. 

• The rate obtained from a stationary state of thermal equilibrium is then assumed to be valid in each fluid 
cell (of spatial size larger than the mean free path) which is taken as the local rest frame. The invariant rate 
(independent of time) is then written in terms of the proper time and fluid rapidity by performing a Lorentz 
boost and assuming that the temperature is a function of the proper time. The resulting rate is then integrated 
during the space-time history of the plasma in combination with a hydrodynamic description of the expansion. 
During the hadronization transition the yield is obtained from a Maxwell construction of the coexistence region 
(under the assumption of a first order transition) . The lever rule is invoked to obtain the photon yield from the 
mixed phase. Thus despite the fact that the rate has been obtained by taking the initial and final times to =Foo 
it is used to extract the photon yield during a finite lifetime, and even during phase coexistence | id ITTl Il2| . 

Caveats: 

The main reason that we delve on the specific steps of the usual computation and on the detailed analysis of the 
main assumptions is to emphasize the inconsistencies in applying this approach to an expanding QGP of finite lifetime. 
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• Hydrodynamic evolution is an initial value vroblem\21i |22j , namely, the state of the system is specified at an 
initial (proper) time surface to be of local thermodynamic equilibrium at a given initial temperature, and the 
hydrodynamic equations are evolved in time to either the hadronization or freeze-out surfaces if the equation 
of state is available for the different stages. The calculation based on the S-matrix approach takes the time 
interval to infinity, extracts a time-independent rate treating the QGP as a stationary state, and inputs this 
rate, assumed to be valid for every cell in the comoving fluid, in the hydrodynamic evolution during a finite 
lifetime. 

• There is also a physical inconsistency in using the S-matrix yield in a hydrodynamic evolution for very large 
photon energy. A hydrodynamic description, which is based on local thermodynamic equilibrium, is valid on 
spatial scales larger than the mean free path for parton-parton collisions in the plasma A ~ 0.5 fm. Thus photon 
momenta k > 2 — 3 Gev probe distances shorter than the mean free path, and most likely the contribution to 
the direct photon spectrum for transverse momenta larger than about 2 — 3 Gev cannot be reliably extracted 
from a S-matrix calculation. Recent measurements of elliptic flow at RHICjT^ suggest that hydrodynamics is 
a reliable description up to pt ~ 2 Gev/c but the data for V2(pt) show larg e departures from hydrodynamics 
(including pQCD)for px > 2 Gcv/c for charged particles (minimum bias) 23J. 

The low energy region of the photon spectrum is dominated by pion decay and bremsstrahlung in the hadronic 
phase and after freeze-out. Thus, the (transverse) momentum interval in which direct photons could be reliable 
experimental probes of a thermalized QGP is 0.1 Gev < k < 3 Gev. 

• Yet another caveat is that despite the fact that the time interval is taken to infinity, namely much larger than the 
photon thermalization time scale, photons are assumed not to thermalize and to leave the medium. The buildup 
of the photon population is neglected under the assumption that the mean free path of the photons is larger than 
the size of the plasma and the photons escape without rescattering. This assumption also neglects the prompt 
photons produced during the pre-equilibrium stage. Indeed, Srivastava and Geiger [24| have studied direct 
photons from a pre-equilibrium stage via a parton cascade model that includes pQCD parton cross sections and 
electromagnetic branching processes. The usual computation of the prompt photon yield during the stage of a 
thermalized QGP assumes that these photons have left the system and the computation is therefore carried out 
to lowest order in a em with an initial photon vacuum state, namely also ignoring the virtual cloud of photons 
that dress the charged particles in the plasma. As emphasized above, in taking the final time f/ to infinity in 
the S-matrix element the assumption is that the thermalized state is stationary, while in neglecting the buildup 
of the population the assumption is that the photons leave the system without rescattering and the photon 
population never builds up. These assumptions lead to considering photon production only to the lowest order 
in a em , since the buildup of the photon population will necessarily imply higher order corrections. Although 
these main assumptions are seldom spelled out in detail, they underlie all S-matrix calculations of the photon 
production from a thermalized QGP. 

• All calculations of the rate based on the S-matrix approach, obtain the imaginary part of the photon polarization 
to lowest order in a em and in a perturbative expansion in terms of a s (including leading logarithmic terms). 
This expansion assumes that a s is small but the coupling depends on the temperature scale. While it could be 
argued that at the initial temperatures expected to be achieved at RHIC and LHC a s may be small, clearly 
the perturbative expansion breaks down in the expanding scenario, when the temperature becomes near the 
critical for hadronization T c ~ Aqcd ~ 160 Mev. Thus, the regime of validity of all S-matrix calculations of 
the yield in a perturbative expansion in a s is actually limited by the lifetime of the QGP. Hence a perturbative 
evaluation of the rate must be understood to be valid on a time scale of the order of or shorter than the actual 
lifetime of the QGP, which then casts further doubts on the infinite time limit. 

As stated in the introduction, however, the QGP produced in ultrarelativistic heavy ion collisions is intrinsically a 
transient and nonequilibrium state. Since the spectrum of direct photons is deemed to be a clean experimental probe of 
the early stages of evolution of a QGP, it is therefore of phenomenological importance to study nonequilibrium effects 
on direct photon production from an expanding QGP with a finite lifetime with the goal of establishing potential 
experimental signatures. 

The current understanding of the QGP formation, equilibration, and subsequent evolution through the quark- 
hadron (and chiral) phase transitions is summarized as follows. A pre-equilibrium stage dominated by parton-parton 
interactions and strong colored fields which gives rise to quark and gluon production on time scales < 1 fm/c |l5j . 
The produced quarks and gluons thermalize via elastic collisions on time scales ~ 1 fm/c. Hydrodynamics is probably 
the most frequently used model to describe the evolution of the next stage when quarks and gluons are in local 
thermal equilibrium (although perhaps not in chemical equilibrium) pH |22|. The hydrodynamical picture assumes 
local thermal equilibrium (LTE), a fluid form of the energy- momentum tensor and the existence of an equation of 



7 



state for the QGP. The subsequent evolution of the QGP is uniquely determined by the hydrodynamical equations, 
which are formulated as an initial value problem with the initial conditions specified at the moment when the QGP 
reaches LTE, i.e., at an initial time U ~ 1 fm/c. The (adiabatic) expansion and cooling of the QGP is then followed 
to the transition temperature at which the equation of state is matched to that describing the mixed and hadronic 
phases E3, |H El ■ 

Our main observation is that the usual computations based on S-matrix theory extract a time independent rate 
after taking the infinite time interval, which is then used in a calculation of the photon yield during a finite time 
hydrodynamic evolution. There is a conceptual inconsistency in this approach, which merits a detailed study based 
on the real time evolution of the photon distribution, which we undertake below. 



III. REAL TIME APPROACH 



A. Gauge Invariance 

Before we focus on the calculation of the photon yield in real time, we address the issue of abelian gauge invariance 
to highlight that the results of the real time approach are fully gauge invariant. Since the relevant interaction is 
electromagnetic, we focus our discussion on the abelian gauge coupling. Let us consider the following Lagrangian 
density for one massless fermion species 



£ = TP(ip-e4)iP--F^F» v , 



(III.1) 



where the zero-temperature mass of the fermion m has been neglected since we consider the high temperature limit 
T m. We begin by casting our study directly in a manifestly gauge invariant form. In the abelian case it is 
straightforward to reduce the Hilbert space to the gauge invariant states and to define gauge invariant fields. This is 
best achieved within the canonical Hamiltonian formulation in terms of primary and secondary class constraints. In 
the abelian case there are two first class constraints: 



0, V 



-e ip^ip 



(HL2) 



where ttq and tv = — E are the canonical momenta conjugate to A and A, respectively. Physical states are those 
which are simultaneously annihilated by the first class constraints and physical operators commute with the first 
class constraints. Writing the gauge field in terms of transverse and longitudinal components as A = Al + At with 
V x Al = ; V ■ At = and defining 



(III.3) 



where G(x — y) the Coulomb Green's function satisfying V^G(x — y) = £ 3 (x — y), after some algebra using the 
canonical commutation relations one finds that At(x) and ^'(x) are gauge invariant field operators. 

The Hamiltonian can now be written solely in terms of these gauge invariant operators and when acting on gauge 
invariant states the resulting Hamiltonian is equivalent to that obtained in Coulomb gauge. However we emphasize 
that we have not fixed any gauge, this treatment, originally introduced by Dirac is manifestly gauge invariant. The 
instantaneous Coulomb interaction can be traded for a gauge invariant Lagrange multiplier field which we call A , 
leading to the following Lagrangian density 



: ¥(* - e 7 A° + ej ■ A T )* + )- [(d M A T ) 2 + (VA ) 2 ] 



(III.4) 



We emphasize that A should not be confused with the temporal gauge field component. 

In the gauge invariant sector of the Hilbert space, namely between states annihilated by the first class constraints, 
and generalizing to Nf flavors of quarks, the Hamiltonian is given by, 



H 



d 3 ;> 



B 2 



&(-ia- V)\P + e J - A T 



(III.5) 



where 



(III.6) 
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is the gauge invariant current and H cou i is the abelian Coulomb interaction which will be irrelevant for our consid- 
erations. While this Hamiltonian is equivalent to that obtained in Coulomb gauge we emphasize that we have not 
imposed any gauge fixing, the Dirac procedure is manifestly gauge invariant. In particular the interaction part of the 
Hamiltonian that will be relevant for the discussion of photon production to lowest order, namely 



Hi = e J df x J ■ A T (III.7) 

is manifestly gauge invariant. Furthermore, states constructed out of the non-interacting Fock vacuum by combinations 
of the gauge invariant operators At, Et are obviously gauge invariant. 

This discussion makes explicit the gauge invariance of the real time formulation. 

Including now the non-abelian color interaction between quarks and gluons the total Hamiltonian is given by 

H = Hqcd[V] + J d 3 x i (E T + B 2 ) + e J d 3 x J • A T + H coui , (IH.8) 

where Hqcd is the QCD Hamiltonian in absence of electromagnetism but in terms of the gauge invariant (under 
abelian gauge transformation) fermion field W and the subscript T refers to transverse components. We have extended 
the fermion content to Nf flavors and the charge of each flavor species in units of the electron charge is included in 
the corresponding current. 

Expanding the gauge invariant field At in terms of creation and annihilation operators in a volume V 

A T (x) = J2 - J= Uk,A e lk ' x + 4 a e-H , (HI.9) 
k,A ^ 2Vk 1 J 

with ek,A the usual transverse polarization vectors, the gauge invariant photon number operator for polarization A is 
given by 

™k,A = a k.A a k,A • (III. 10) 



B. Time evolution 



Time evolution in quantum mechanics and quantum field theory is an initial value problem. Given the total 
Hamiltonian H the time evolution of a density matrix is completely determined by specifying the density matrix at 
an initial time to. Once the initial density matrix is specified at the initial time, its time evolution is completely 
determined by the unitary time evolution operator e - 4 ^ f (*-*o) > j n the case under consideration the Hamiltonian H is 
time independent and gauge invariant and given by eq. (|III.8|) . The ensuing real time dynamics is therefore completely 
specified by prescribing the initial density matrix at the initial time to. 

As discussed above the S-matrix calculation implicitly assumes that the initial density matrix is that of thermal 
equilibrium for quarks and gluons at an initial time to — — oo and explicitly assumes that initially there are no photons. 
With this choice of initial conditions, since the initial state has been prepared at to — — °o there is no memory of the 
initial state, in agreement with a stationary state in thermodynamic equilibrium. 

Our goal is to relax this assumption of a thermal stationary state and study the consequences of a true real time 
evolution as befits the physical problem of a thermalized QGP emerging about 1 fm/c after a nucleus-nucleus collision 
and evolving during a finite lifetime of about 10-20 fm/c towards a hadronization (and confinement) transition. 

The real time evolution of the density matrix after the initial time to is given by 

p(t) = e - iH{t - to) p(t ) e lH{t ' to) . (mil) 

The total number of photons of momentum k per unit volume at a given time t, namely the photon yield is given by 
(assuming translational invariance) 

^ 3 -dWk = ^ Tr[m ^ x] ' (IIL12) 

A 

with n^A — a A- I n order to compute the direct photon yield to lowest order in the electromagnetic coupling, 
it is convenient to write the total Hamiltonian given by ea. (|III.8ll as 

H = H + Hi, (III. 13) 

Ho = H QCD m + J d3x l (Et + B 2 ), Hi = J d 3 x J • A T , 
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where the current J is given by eq. (|III.6p and we have neglected the Coulomb term since it will not contribute to 
the direct photon yield to order a em . 
The time evolution operator is given by, 

e -iH(t-t ) = e -iH t u ^ ^ e iH t ( (111.14) 
where U(t, to) the unitary time evolution operator in the interaction picture of Hq given by 

U(t, t ) = 1 - i f dt' Hi(t') + 0{e 2 ) ; #/(t) = e lHat ff, e - lHot , (111.15) 
Jt 

where we have only considered the lowest order in the electromagnetic coupling. 

It is then convenient to pass to the interaction picture of Ho namely the full QCD Hamiltonian and free electro- 
magnetism by defining the initial density matrix in the interaction picture of Ho as 

p lp (to) = e lH ° t0 p{to) e- m ^ . (III. 16) 

For the case of interest the initial density matrix describes a quark gluon plasma in equilibrium under the strong 
interactions, it must therefore commute with Ho in which case Pi P (to) = p(to). At any given time t 

p{t) = e- lHot U(t, t ) p ip {to) U-% t ) e iHot . (111.17) 

Since the number operator a£ A a k ,\ commutes with Hq we find for the direct photon yield at time t the following 
exact expression 

(2tt) 3 = £ Tr [U(t, to) p lp (t ) U-% t ) n k)A ] . (111.18) 

In the interaction picture of Ho the time evolution of J and At is given by 

3(x,t') = e tH ^ CDt ' 3{x) e - iHc * CDt ' , 
£k,A 

k.A 



A T (2 , = ]T [a k>A e* k ' x + 4, A e^ * e M '\ , (111.19) 



where e k ,A ar e the transverse polarization vectors. 

Given an initial density matrix, the photon yield can be calculated by inserting a complete set of eigenstates of H$ 
and computing the required matrix elements. 

We note that the real time expression for the photon yield eq. lIII.lSp is exact in terms of the full time independent 
Hamiltonian and is gauge invariant provided that the density matrix is constructed with physical states. Evolution 
of quantum states or density matrices is an initial value problem, once the density matrix has been specified at an 
initial time, its time evolution is completely determined by the Hamiltonian. 

We carry out this program below to lowest order in the electromagnetic coupling for cases that are relevant to the 
description of direct photons from a QGP. 



IV. PHOTON YIELD 



We compute here the photon yield for a QGP with finite lifetime in the real time approach. 



A. Thermalized QGP, no initial photons: 



We begin the study of the real time dynamics of direct photon production by considering that at an initial time 
to ~ 1 fm/c the QGP is in thermal equilibrium (under the strong interactions) and there are no initial photons. Such 
initial density matrix is compatible with the usual assumption on the initial state invoked in the S-matrix calculation 
described in section [H] but makes explicit that such state describes an initial value problem from a thcrmalization 
time to taken to be of the order of 1 fm/c. 

p(t ) = J2^ En " K>NI ® IO7XO7I ; H QCD \n q ) = E nq \n q ) , (IV. 1) 

n q 
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and |0 7 ) is the photon vacuum, annihilated by the gauge invariant operators ak,A- The incoming nuclei coasting along 
the light cone are exact eigenstates of the full Hamiltonian H. Hence, neglecting bremsstrahlung off the Coulomb 
field of the nuclei there is no photon emission prior to the collision. Thus, the state |0 7 ) is the 'in' state up to the 
time of the collision, which is annihilated by the 'in' operator a k \. Therefore the initial density matrix eq. I|IV.1|) is 
consistent with neglecting the photons produced between the collision and the onset of thcrmalization, namely during 
the pre-equilibrium stage. 

We will study in detail alternative and more general initial states that include photons and correlations in a later 
section (see section Ivl below) . 

Obviously the initial density matrix given by ea. (|iV.l| l is gauge invariant and commutes with Hq in ea. l|III.14|) . it 
follows that p(t Q ) = Pi p (t ) and it describes a thermal ensemble in equilibrium under the strong interactions with no 
photons. This assumption is compatible with all the calculations of direct photon production from an equilibrated 
QGP available in the literature. Studying the time evolution of this initial state allows us to address the dynamics of 
the formation of the virtual photon cloud and to highlight the inherent difficulty in separating the observable photons 
from those in the virtual cloud in the plasma during a finite lifetime. 

Defining as |ti 9 ;to 7 ) = \n q ) ® |ni 7 ) the eigenstates of Hq with m 1 photons (we do not specify the wavevector and 
polarization to avoid cluttering of notation) we can compute the matrix elements in ea. (|III.18|) by inserting a complete 
set of these eigenstates. To lowest order in the electromagnetic coupling the set of intermediate states that contribute 
to the photon number (hk,\) contain only one photon of momentum k and polarization A. 

After a straightforward calculation we find, to lowest order in the electromagnetic coupling 



dN(t) _ e 2 



-ik(t 2 — ti) ik-(x 2 -xi) 



Writing, 



Y e~ 0En * (n 9 |ek,A • J(x 2 ,t 2 )\m q ) (m 9 |e kiA • J(x[,ti)\n q ) . (IV.2) 

n q ,m q 

3{x, t) = e- lp - 3 e %H » CDt J(0, 0) e ^ H » CDt e ips , (IV.3) 



with P the momentum operator and choosing the eigenstates \n q ) to be simultaneous eigenstates of P with P\n q ) = 
Pn q \n q ), we can write 

]T e- 0E ^ (n q \e k ,x-J(x-2,t2)\m q ) (m q \e k , x ■ J(afi,ti)K) = f d 3 p dcu e -tf-(a.-*0 e M* a -*0 4^ ^ a^(p,u>) , 

n q ,m q 

(IV.4) 

where 

a>-{p,u) = ]T e-P E "« (n g |Ji(0,0)|m g ) (m q \ J,(0, 0)\n q ) 6 3 (p~ (p nq - p mq )) S(to - (E riq - E mq )) . (IV.5) 

n q ,m q 

Carrying out the integrals in ea. l|IV.2|) and summing over the polarizations, we finally find, 

dN(t) e 2 f, . 1 - cos[(w - k)(t - 1 )] ,„ T 

where 

(ft (IV.7) 

is the transverse projector. 
In appendix lAl we show that 

(2i) 3 e 2 <7>(^)^Imn tJ (^) ; n(u) = -^—j, (IV.8) 

where ImIl i j(k 1 Lu) is the imaginary part of the retarded photon polarization tensor. However, the photon production 
yield is calculated to lowest order in a em and in principle to all orders in a s . 
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Introducing the transverse photon polarization 

n T (k,u J )~v^(k)n i] (k,Lj) , (iv.9) 

we finally find the real time expression for the photon yield to be given by 

dN(t) _ 1 f +oc doj Imn T (k,u) l-cos[(w-fc)(t-*o)l ay 10) 



d 3 xd 3 k (27r) 3 k J tt eT — 1 (w — fc) 2 

where the photon polarization is obtained in lowest order in a em and in principle to all orders in a s . This expression 
coincides with the lowest order result obtained from a kinetic description in ref.|20j|. 
If the long time limit is taken and the following identity is used (see poH^ 1 ) 

«• ^ i-co^-^-M ,-„,^ <t _ + p ^ „ i + fi\ 

-oo (uj-ky {uj-ky \t-hj 

then if ImiTT(fc,w — k) ^ 0, we obtain the long time limit of the total yield, 

dN(t) _ i f imn T (fc,c = fc) (f ^ | pd. imn T (^) p i | 2) 



(IV. 13) 



d 3 a:#/c (27r) 3 fc 1 e #_i v uy J.^ tt e* - 1 (u - fc) 

^From the result (|IV.12(I above, we obtain the long time limit of the invariant rate 

dN(t) t-to^+oo 1 ImII T (fc,a; = fc) 
d 4 xd 3 k ~ (2tt) 3 e * _ l 

This is the same as the S- matrix result given by ea. (|II.8|l . the expression (IIV.12II only involves the transverse part of 
the photon polarization, a consequence of the manifestly gauge invariant treatment. Thus, we highlight that in the 
infinite time limit our result for the invariant rate coincides with the usual one. 

However, taking the time to infinity introduces all the caveats that were discussed in section [H] above. Thus, in 
order to avoid these caveats, the initial time to must be interpreted as the time at which quarks and gluons thermalize. 
The expression (|IV.10f) determines the total photon yield (per unit volume) at a given time t, thus the total direct 
photon yield at the hadronization time is obtained by setting t = thad ~ 10 — 20 fm/c for RHIC and LHC respectively. 

The asymptotic long time expression (|IV.12|I . features two different terms. The first term, which grows linearly in 
time leads to the rate which is time independent and is associated with the photons that are produced per unit time 
per unit volume. The second, time independent term can be interpreted as the total number of photons in the virtual 
cloud dressing the quarks in the medium. 

This manifest separation between the time independent term that describes the photon cloud associated with the 
charged particles, and the photons produced at constant rate (namely the term in the yield that grows linearly with 
t — to) only emerges in the long time limit. 

For any finite time interval there are contributions to the photon yield from the whole range of u> ^ fc. The long 
time asymptotic behavior of the total photon yield (|IV.10p is determined by the behavior of ImiT(w ~ fc). In particular 
if ImII(fc,c<j = fc) = there is no linear time dependence asymptotically, however if ImII(w ~ fc) tx (u> — fc) as is the 
case for the one loop contribution to the polarization [of 0(a ern a s ), see below] then using the formula [20l l27j 

dy^L [1 - cos(yfci)] *~'= DO p(0) m(fci) + 0(1) , (IV.14) 

o y 

for p(oo) = 0. We thus conclude that the photon yield will grow logarithmically in time in this case. 

The long time behavior of the yield is therefore determined by the behavior of the ImIT(fc, u>) in the region ix> ~ fc. 
Generally, the imaginary part of the photon polarization for oj ~ fc behaves as 

Imn T (fc, uj) "= fc Imn T (fc, lu = fc) + ImiT^fc, lu = fc) (u - fc) + O [(cj - fc) 2 ] , (IV.15) 
then we findpoll^ the long time behavior of the yield to be generally given by 
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where the prime in the second term stands for a derivative with respect to the frequency. Obviously at very long times 
the term with the linear time dependence will ultimately dominate. However, for a finite time interval (t — £ ) 7 it is 
possible that a lower order logarithmic term may give a comparable contribution to or even be larger than a higher 
order term linear in time. Furthermore, it is clear from the above analysis that the contributions to the imaginary part 
of the photon polarization for uj ^ k will actually contribute to the yield during a finite time interval. However, these 
contributions are absent in the S-matrix calculation which only extracts the linear time dependence in the asymptotic 
long time limit given by ea. HIV. 130 . 

The main point of this discussion is to highlight the following important issues: 

• The real time calculation reproduces the result of the S-matrix approach in the asymptotic long time limit. 
Therefore the usual S-matrix result is contained in the real time approach, which provides a detailed description 
of the process of photon production during a finite lifetime. 

• During a finite time interval the region lu =/= k of the imaginary part of the photon polarization contributes. 
The contributions from different regions of the spectral density have different time dependence. During a finite 
time interval contributions that are subleading in the asymptotic long time limit can be of the same order as 
the term that becomes asymptotically linear in time and which defines the rate. In a finite time interval the 
contributions from the different regions in the frequency integral cannot be separated. Thus, in order to reliably 
understand the time dependence of the yield during a finite time interval, one must find the imaginary part of 
the photon polarization in the full range of frequency, not just at lu = k which only determines the asymptotic 
long time behavior. 

• This discussion makes manifestly clear that the photon yield obtained from the S-matrix calculation of the rate, 
namely extracting the linear time dependence in the asymptotic long time limit, ignores all other contributions 
which grow slower in time but that do contribute to the yield for a finite time interval. 



B. Photon production in the hard thermal loop approximation 



To make the above statement more quantitative and to begin our study of the real time description of photon 
production within a specific example highlighting the conclusions above, we beg in by considering the imaginary part 
of the photon polarization tensor in the hard thermal loop approximation[28Tl29| . This approximation yields the 
leading result for the polarization tensor for soft photons, namely t«r and arises solely from in-medium processes. 
While we will obtain the full one-loop contribution to the imaginary part of the photon self energy below, the HTL 
limit only features medium dependent contributions. Thus the HTL limit allows us to address dynamical issues solely 
associated with in-medium processes. 

For two flavor of quarks (up and down, with three colors), the imaginary part of the transverse photon polarization 
in the hard thermal loop approximation is given bv|20| 



ImH 



HTL 



407r 2 a em T 2 LU ( LU 2 . 2 , 



(IV. 17) 



Obviously this contribution to the imaginary part vanishes linearly as u> — > k. However, while ea. (|IV.12|) would lead 
to a constant, time independent yield and therefore to a vanishing invariant rate since lmHij(uj = k) = 0, the correct 
asymptotic long time limit of the yield follows from ea. l|IV.16|) and is given bvp(i| 



dN 



HTL 



(*) t-t 



++oo 5 ot em T 



d 3 xd 3 k 



18tt 2 k 2 



n(k) \n[k{t - t )] + — I duj^-V 



2k 



— k 



k 



e T 



1 



O 



t - tr 



(IV. 18) 



This expression would lead to an invariant rate that vanishes as 0(l/t) for t — > +oo. However, for any finite time 
there is a non-vanishing contribution to the yield. 

The second, time independent term in the eq. (|1V.18|) can be identified with the virtual photon cloud in the thermal 
bath. This contribution is medium dependent and can only be identified in the long time limit, the separation between 
the time independent and the time dependent contributions is meaningful only in the long time limit. For any finite 
lifetime there is no unambiguous separation between the different contributions, and only the full photon number is 
meaningful. 

A comparison with the photon equilibrium spectral density in the HTL approximation 29] reveals that the integrand 
of the time independent part in eq. ljIV.18JI is simply related to the HTL wave function renormalization and the real 
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part of the photon polarization. Namely, the virtual cloud is actually revealing the dynamics of formation of a plasmon 
quasiparticle in the medium. 

The result ea. l)iV.18|l is valid in the asymptotic long time limit, for any arbitrary finite time we must use the full 
expression given by ea. (|IV.10|l . In particular, in the hard thermal loop limit we find that the yield is given by 



dN 



HTL 



(t) 



d 3 xd 3 k 



kr, — 

' T 



5 a em T 

36tt 2 k 2 
1 

dx n 

-l 



T 

k~ 
x — 
T 



k(t -to); 



T 



x{l-x z ) 



2 1 — cos [{x — l)kr] 



(x-iy 



(IV.19) 



The function T [kr; =1 is displayed in fig. 2] for several values of the ratio k/T. This figure displays the logarithmic 
growth determined by the asymptotic behavior (|IV.18(I at long times. 




FIG. 1: The function T [kr; |] vs. fcr for £ = 0.1; 0.3; 0.5 respectively 



C. Full one loop polarization: 

As mentioned in the introduction, we will study the dynamics of photon production during a finite lifetime by 
focusing on the lowest order contribution to the photon polarization. This is a quark one loop diagram and is of order 
a em . The imaginary part of this diagram vanishes at u) = k, hence this lowest order diagram does not contribute to 
the usual rate obtained from the long time limit as discussed in detail above. 

The goal of this study is to understand the photon production during the finite lifetime of the QGP from processes 
whose contribution is subleading in the asymptotic long-time limit. Of course diagrams which are higher order in a s 
will also give contributions to the photon production from the region lo ^ k, but focusing on the lowest order diagram 
we will be able to extract important aspects of the dynamics that are missed by the S-matrix calculation and that 
could be experimentally relevant. 

A comparison between real time yield obtained from the one loop contribution to the photon polarization and the 
S-matrix yield in a wide range of photon energy, in particular for large photon momentum, requires the full expression 
for the photon polarization to one loop order. 

A lengthy but straightforward calculation gives the following result 

Imn T ( W , k) = 7T (W, k) + 7T LD (iO, k) + TT 2P (u>, k) (IV.20) 

with 7r ; 7r 2 p and ttlo the zero temperature, two fermion thermal cut and Landau damping contributions, respectively, 
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given by the following expressions 
10 



7T (oj, k) 
7T 2 p(w, k) 



n LD (uo, k) 



9 

10 



ctem (uo — k )Q(uo — k ) sign(cj) 



U3 



In 



1 + e" 



T V 1 + e 



2T 



E(- ir +1 



I —mW- g— mW+ 



Tin' 



(e- mW - 



+ e 



12 T2^_ W 

3 ttem ^ fc 2 

— mW- g — mW + ^ 



Am 

T 



6(w 2 - fc 2 )sign(w) 



1 + e 



1 + e-^+ 



2T 



m 3 Tm 5 



9(fc 2 - uo 2 ) sign(w) , w± 



±k 



2T 



(rv.21) 



(IV.22) 



(IV.23) 



The first two terms 7To,7T2p arise from the process qq — *■ 7 and 7T£,£) from in-medium bremsstrahlung q — > jq. The 
long- wavelength limit k <C T is dominated by 7T££> and simplifies to the HTL expression eq. I|1V.17|) . 



D. Dynamics of the virtual photon cloud 

As discussed above, the asymptotic long time limit is determined by the behavior of Imfly (A;, uo) for uo ~ k, therefore 
it is convenient to separate the contribution from the positive and negative frequency regions in the integral in eq. 
(|IV.10fl . Using the properties ImIlV(fc, —uo) = — LxuT/j^fc, w) and n(— u) = —[1 + n(u>)] we write 

dN(t) = dN(+Ht) dN(-\t) 

d 3 xd 3 k d 3 xd 3 k d 3 xd 3 k 1 ■ ' 

in terms of the positive (+) and negative (— ) frequency contributions respectively, given by 

dN^jt) 1 f°° duo l-cos{(uo-k)(t-t )} 

dN ( -\t) 1 f°° duo , , %1 „ x 1 - cos[(w + ft)(i - i )l 

«f = (2^ /„ v t 1 + ^ imn ^ fc ' — (i +fc ) 2 (iv - 26) 

For further analysis, it is convenient to add together the terms that feature the Bose-Einstein distribution function, 
thus we write instead, 



dN dN^ dN<W 



d 3 xd 3 k d 3 xd 3 k d 3 xd 3 k ' 



(IV.27) 



with 



WW® _ 1 f-du (l-co S [(uo-k)(t-t )} + l-cos[(, + AQ(t-M]. j(i ^ 28) 



d 3 xd 3 k (2ir) 3 kJ ir v ' v ' '{ {uo - k) 2 (10 + k) 

dN(v) ^ - 1 r^^^-^f+m-to)], (rv . 29) 



d 3 a;<i 3 /c (2tt) 3 k J ir v ' ' (uo + k) 2 

Both terms are positive, however, while the frequency integral in d ^p xd 3^ is finite because of the Bose-Einstein 

distribution function, the frequency integral in d ^i x ^ features divergences associated with the virtual photon cloud, 
which can be seen as follows. 

The contribution d ^i xd -})l > does not feature resonant denominators, therefore it remains positive and does not grow 
in time at long time. The oscillatory terms average out on a short time interval, as shown explicitly in fig. [3 which 
displays the negative frequency contribution integrated up to a frequency cutoff uo c = 100 Gev. 

Thus, after a very short transient time the average of d % xd -if) is obtained by neglecting the oscillatory cosine term 
in eq. (H\T29l) . 
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FIG. 2: The contribution d % x J£ given by eq. 11V.29I integrated up to w c = 100 Gev for k = 3 Gev/c; T = 0.3 Gev . 



The remaining time independent term in eq. l|IV.29p features two distinct contributions: 

i) the zero temperature part of Imlly, namely ttq given by ljIV.21J) yields the distribution of photons in the virtual 
photon cloud of the vacuum given by 

i r 



d 3 xd 3 k k J k (uj + k) 2 8vr 4 

d 3 xd 3 k 



This contribution can be extracted by taking T — and to — > — oo in ^jrX , it is linearly divergent and it clearly 



must be subtracted since it is not observable. 

ii) The finite temperature contributions from 1T2P and ttld for k 3> T are dominated by the region u> ~ k in the 
frequency integral in eq. (jIV.29|) which is finite. A lengthy but straightforward calculation from eas. (|IV.22|) and 
(IIV.23|I yields the following result for the frequency integral of both these terms for k T 



d 3 xd 3 k 



dN vv _ 10cwC(3) T 3 , /J_ . (1V , S]) 



d 3 xd 3 k 32 7T 4 fc 3 \k 5 



where 'av' refers to the time average of the oscillatory term and £(3) = 1.202057. . .. The result (|IV.31|) originates 
in the region lu ~ — k both in the Landau damping as well as two particle contribution. This region features a 
contribution that is not exponentially suppressed in k for k^> T . 

It is clear that the integral over the momenta k yields a logarithmically divergent number of photons in the medium 
and a linearly divergent energy integral. We identify this temperature dependent term as describing the virtual photon 
cloud in the medium, which again must be subtracted since it is unobservable. 

The divergence associated with the virtual cloud of the vacuum given by eq. Ijl V .301) as well as the divergences in 

the photon number and energy stemming from the ea. i|iV.31|) must both be subtracted from d % xd 3)} ■ 

An alternative and illuminating interpretation of the virtual cloud of the vacuum emerges by noticing that eq. I|IV.30|) 
is related with the vacuum wave function renormalization Z given by 



Re 



dnT(oj,k) 



duj 2 



1 f + °° , lmU T (uj',k;T = 0) 



j=k 



^' , ' n 2 . (IV.32) 



2nkJ_ OQ (k-Lu') 



where we used the dispersion relations eas. (|A.6(l - HA.7|) . We see that ea. HIV. 30(1 is just the zero temperature wave 
function renormalization to orden e 2 . The vacuum virtual cloud dresses the bare particles into the physical 'in' or 
'out' states, the wave function renormalization is simply the overlap of these. 
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The asymptotic reduction formula (LSZ formulation) requires that 



ak,A ; 4,A ~> Z "a k:A;out ; Z 2 4,A;out (IV.33) 

Thus multiplying the number operator by Z~ x , cancels the vacuum contribution to 0(e 2 ) thus justifying the 
subtraction of the vacuum term. The asymptotic reduction (LSZ) formalism requires the zero temperature wave 
function renormalization since the in and out states are the states created from the physical vacuum by the in and 
out operators. 

In the medium the extra in-medium contribution to the virtual cloud dresses the physical particle into a quasi- 
particle, in this case a plasmonj^l these are not asymptotic states. Thus the time dependent terms which are 
associated with the virtual cloud at asymptotically long times are actually describing the dynamics of formation of 
the quasiparticle in the medium. 

While subtracting the vacuum term corresponds to multiplying by the inverse of the zero temperature wave function 
renormalization according to the LSZ reduction formula and the propagation of a physical particle in the out state, 
a similar interpretation for the medium contribution is not available. The virtual cloud of the medium dresses the 
physical particle into a quasiparticle as it evolves in the medium but is not an asymptotic state. The subtraction of 
the in medium virtual cloud cannot be justified on the basis of asymptotic theory. 

Furthermore, there is no unambiguous manner to subtract the divergent terms at all times, since the oscillatory 
terms are finite, and subtracting solely the time independent terms leaves an expression that becomes negative, 
features the divergences described above at the initial time t — to but is finite for any t ^ to. 

Hence, in this section we proceed to subtract completely the contribution d % xd Ai} an d use the following definition 
of the subtracted yield for the analysis that follows. 



dN{t) = dN^ T \t) _ 1 f°° duj 
d 3 xd 3 k = d 3 xd 3 k ~ (2tt) 3 k 



f duj 

/ — ntoj) ImlLr (k, to) 
Jo f 



1 — cos(cj — k)(t — t ) 1 — cos(ix) + k)(t — to) 



(w - k) 2 (w + kf 



(IV.34) 

This expression is finite and positive at all time and the large lo, k regions are exponentially suppressed. 

We emphasize that completely subtracting d % xd 3^ also neglects the positive and finite time dependent contributions 
from this term but which cannot be unambiguously separated from the divergent terms during a finite time interval. 

At asymptotically long time the definition of the yield l|IV.34|l features a term that is constant in time and terms 
that grow either linearly or logarithmically (or both) in time. The time independent and finite contribution which is 
typically associated with the virtual photon cloud can be separated unambiguously from the photons produced with 
a constant rate only in the asymptotic long time limit. However, since the QGP has a finite lifetime, only the virtual 
cloud of the vacuum and the medium contribution that leads to a divergent number of photons and energy can be 
unambiguously associated with the virtual photon cloud of the medium. Furthermore, these (divergent) contributions 
are associated with very fast oscillations and become constant in time on a very short time scale (<SC 1 fm/c) as clearly 
shown in fig. |21 

To be sure the photons produced by the plasma are detected far away from the collision region and (practically) at 
infinite time. However, these photons had been produced in the plasma during a (much shorter) time scale to < t < tf 
with tf being the hadronization time. Clearly the vacuum part of the virtual cloud can be recognized and subtracted 
unambiguously, it is given by eq. (|IV.30|) . However, identifying the contribution to the virtual photon cloud in the 
medium can only be achieved unambiguously if the formation time of the virtual cloud is much shorter than the 
lifetime of the QGP and the number of photons in the cloud diverges. As mentioned above the time dependent 
terms which asymptotically are associated with the virtual cloud are describing the dynamics of formation of the 
quasiparticle in the medium. 

In section [V] below we will introduce and implement a method that allows to separate the divergent contributions 
to the virtual cloud which are responsible for the rapid oscillations, in an effective manner. 

The contribution ^p-grjr evolves on longer time scales and contains all of the potentially secular terms, those 
that grow linearly, logarithmically etc. Furthermore, since the large lo regions are exponentially suppressed, the 
integrals are dominated by the region lo ~ k and lo ~ 0. The region to ~ k leads to secular terms from the resonant 
denominators, the non-resonant terms feature oscillations on the time scale ~ 1/fc, hence their time dependence is 
relevant during the lifetime of the QGP specially for long wavelengths. Thus, no subtractions on this contribution 
are warranted. 
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E. Real time vs. S-matrix yields: 



We can now establish a comparison between the yields and spectra predicted by the real time expression eq. I|IV.34|) 
with those obtained from the S-matrix approach (the equilibrium rate). The equilibrium rate in leading logarithmic 
order approximation in the strong coupling a s is given in ref.@, and for two flavors (up and down) of quarks (and 
three colors) becomes 
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, (IV.35) 



where rif(k) is the Fermi distribution function. This fit seems to be very accurate in the region of momenta 0.2 < 
/c/T<50|3. 

We will also use the lattice parametrization [20j for the temperature dependence of the strong coupling a s (T) given 
by 



<**(?) = 



Grr 



29 In 17 



T c ~ 0.16 Gev 



(IV.36) 



Although this lattice fit is valid at high temperatures and certainly not near the hadronization phase transition, wc 
will assume its validity in the temperature range relevant for RHIC in order to obtain a numerical estimate of the 
S-matrix yield. We note, however, that at a s (0.3 Gev) ~ 0.24 and the validity of the perturbative expansion is at 
best questionable. 

The yield per unit phase space as a function of time obtained from this rate is given by 
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(IV.37) 



We begin by comparing the yield from the real time evolution with the photon polarization in the hard thermal loop 
approximation, valid for k -C T. 

Since the hard thermal loop (HTL) limit eq. (|IV.17|) is valid for k <C T and the formula for the rate l|IV.35|) is valid 
for 0.2 < k/T < 50 the comparison between the two is reliable for k/T ~ 0.2. 

A comparison of the yields per unit phase space from the time dependent real time expression (|iV.10H given 
by ea. i|fV.19|) in the HTL approximation, and the leading order result of the S-matrix formulation obtained from 
ea. HIV.35|) 8] for k = 0.1 Gev; T = 0.5 Gev is displayed in figureEJ 

It is clear from this figure that during a finite time interval compatible with the expected lifetime of QGP in local 
thermal equilibrium at RHIC, the hard thermal loop contribution is of the same order as the S-matrix result. 

Fig. 21 compares the real time yield ea. ljiV.34() with the full one loop photon polarization given by eqs. (|IV.201 
IIV.23|) to the S-matrix yield for k — 3 Gev/c ; T = 0.3 Gev. The real time yield is dominated by the Landau damping 
contribution and clearly competes with the S-matrix result during the lifetime of the QGP at RHIC. 

Fig. [5] displays the logarithm of the real-time photon yield given by eq. I|IV.28|) vs. k compared to the S-matrix 
result for T — 0.3 Gev ; t — to — 10 fm/c. Cleary both spectra fall off exponentially, but the spectrum from the real 
time yield falls off slower and displays an excess of photons as compared to the equilibrium one for k > 2.2 Gev/c. 
Since eq. (|IV.28|) is dominated by the Landau damping contribution to the imaginary part of the photon polarization, 
these photons originate in bremsstrahlung which is a medium effect. Wc emphasize that we have only considered the 
contribution from eq. I|IV.28() since we have subtracted the full contribution from eq. (|IV.29|) which is not exponentially 
suppressed but falls of as a power law. The subtraction of ea. ljIV.29l) was motivated by the fact that this contribution 
features both vacuum and in-medium divergences associated with the unobservable virtual photon cloud and the finite 
contributions cannot be unambiguously separated during a finite interval of time. We will revisit this point in section 
below where we introduce a formulation that allows to separate the virtual cloud by allowing an initial preparation 
stage. 



F. Energetics 



If the thermalized plasma is emitting photons, the energy radiated away must be drained from the plasma. In this 
section we study the different contributions to the energy radiated away with the escaping photons. The total energy 
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FIG. 3: Comparison between the S-matrix yield SM given by eas. ljIV.37(l - l|IV.35|l and the real time yield (HTL) given by 
ca. llV.19t in the HTL approximation, both for k — 0.1 Gev; T — 0.5 Gev as a function of t to (in ftn/c). jV(fe,£) — d^kcfix ' 
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FIG. 4: Comparison between the S-matrix yield NsM{k,t) given by eas. lllV.35ll - lllV.37t and the real time yield N r t(k,t) given 
by ea. llV.34H for k = 3 Gev/c; T = 0.3 Gev as a function of t — to (in fm/c). The real time yield is dominated by the Landau 
damping contribution given by eq. 11V.23II . N(k,t) = d i k % x ■ 



is conserved, namely 



Tr [p{t)H] = Tr [ e -«*(*-*°) p(t ) e tH ^- to) h] = Tr [p(t Q ) H] , 



(IV.38) 



where H is the total Hamiltonian (|III.13Jl . Passing to the interaction picture of Ho in ea. ljlll.13jl the equation above 
(|IV.38j) becomes 



Tv[p(t ) H]=Tr[p lp (t )U- 1 (t,t ) (H + Hi(t)) U(t,t )] 



(IV.39) 
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FIG. 5: Comparison of the spectra, lniV(ft, t) vs. k between the S-matrix Nsm and real time N r t yield given by ea. l|IV.34(l for 
T = 0.3 Gev and t - to = 10 fm/c. 



where U(t, t ), -ffj(i) and pi P (t ) are given by eqs. I|III.15|) and l|III.16ll respectively and U 1 (t,t ) — U(t ,t). Assuming 
translational invariance in a (large) volume V, the statement of conservation of energy of eq. ljIV.39|) becomes 



A£ QCD (t) + £ 1 (t) + £i(t) = 0, 



(IV.40) 



where the second term above corresponds to the energy per unit volume radiated away by the photons produced in 
the plasma [see eq. <|III.12|) ] and we have introduced the following definitions, 



A£ Q c D {t) = y Tr [p ip (t ) U- l {t,t Q ) H Q cd U{t,t )] - - Tr[p ip (t ) H QCD } , 
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(IV.41) 
(IV.42) 



The individual terms above can be computed to lowest order in a em by expanding the time evolution operator up to 
second order in the interaction. A lengthy but straighforward computation using the initial density matrix given by 
eq. <|1V.1|) and introducing an intermediate set of states leads to the following result 
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(IV.43) 



which manifestly satisfies conservation of energy as in eq. (|IV.40|I . The asymptotic long time limit is determined by 
the region of the spectral density oj ~ k, thus it is clear from the expressions in eq. (|1V.43|I that the interaction energy 
shuts-off in the long time limit and for asymptotically long time the rate of radiative energy loss by the photons is 
balanced by the rate of energy loss of the plasma, which can be interpreted as radiative cooling, namely 



d£~f(t) _ d£ QC D(t) 



dt 



dt 



(IV.44) 



We emphasize that this result is only valid in the long time limit, during a finite time interval there is a contribution 
from the interaction energy which necessarily is present to satisfy energy conservation. 

As discussed in section IIVDI for the virtual photon cloud, the negative frequency contribution features the zero 
temperature divergence associated with the virtual photon cloud of the vacuum, as well as the divergence associated 
with the virtual photon cloud in the medium. The energy in the photon cloud in the vacuum diverges as the fourth 
power of a cutoff, while the energy in the medium contribution of the virtual cloud diverges linearly with a cutoff 
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since the number of photons diverges only logarithmically, as can be seen from eq. (|IV.31f) . Furthermore since the 
virtual cloud builds up in a very short time scale, the negative frequency contribution without the Bose-Einstein 
distribution averages to a time independent constant on a time scale <§; 1 fm/c, as can be gleaned from fig. 
Since this contribution features the divergences associated with the vacuum and in-medium virtual photon clouds, we 
subtract it from the energy, consistently with eq. l|IV.34ll . 
Therefore we now study the following subtracted energies, 
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(IV.45) 



In subtracting the negative frequency contribution without the Bose-Einstein distribution function, we are also ne- 
glecting the finite parts of the negative frequency contribution to £ 7 (i), which as mentioned above cannot be extracted 
unambiguously. 

Fig. El displays the subtracted energy radiated in photons [see eq. (IIV.45I ] as a function of time for T = 0.3 Gev 
as compared to the energy obtained from the S-matrix yield. The real-time energy reveals the logarithmic growth in 
time and it is of the same order as that obtained from the S-matrix yield during the lifetime of the QGP expected at 
RHIC. 




t-t„ (fm/c) 

FIG. 6: Comparison between the energy density radiated in photons as a function of time for T = 0.3 Gev. The real time 
energy radiated in photons E rt (t) displays the logarithmic time dependence. The energy obtained from the S-matrix yield is 
E SM {t). 



Figures[7|display the subtracted contributions eq. (|1V.45|) to AEqcd^) and £i{t), as a function of i for T = 0.3 Gev. 

We confirmed numerically that the main physical mechanism of radiative energy loss is Landau damping by studying 
separately the different contributions to the photon polarization in the energy. 

It is clear from the figure that the interaction term evolves on time scales of order ~ 7 — 8 fm/c, and that the QGP 
cools by photon emission faster than the change in interaction energy, a clear signal that the photons being emitted 
are mainly a consequence of the cooling of the QGP. The numerical analysis reveals that the energy in photons grows 
and that in the QGP diminishes logarithmically, while the interaction energy slowly approaches a constant over the 
time interval of the same order as the lifetime. 
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FIG. 7: Positive frequency contribution to AEqcDit) (left panel) and £i(t) (right panel) vs. t — to for T = 0.3 Gev. 



V. MODELLING THE INITIAL STATE: QUASI- ADIABATIC INTERPOLATING STATES: 

The expression for the total direct photon yield at a given time t given by ea. ljIV.10fl was obtained from an initial 
density matrix corresponding to a QGP in thermal equilibrium and the photon vacuum. The assumption on the 
photon vacuum is in agreement with the assumption of no photons in the initial state in the S-matrix calculation. 
The rational for this choice is that the photons produced from the pre-equilibrium state leave the system without 
buildup of the photon population, and the real time calculation described above makes explicit this choice of initial 
state. 

The study of the dynamics of the photon cloud above indicates that a real time description of photon production- 
a necessary treatment to address the finite lifetime of a QGP- must include the analysis of the virtual photon cloud 
and a consistent and systematic separation between the contributions from the virtual photon cloud and the observed 
photons. 

The S-matrix formulation bypasses any discussion of the dynamics of the photon cloud by initializing at ti — * — oo, 
taking the final time tf — * +oo and extracting only the contributions that lead to a linear time dependence in the 
yield. In this manner, all constant contributions, such as that of the photon cloud, as well as those that grow slower 
than linear in time, are neglected. This can also be summarized with the statement that in the asymptotic long time 
limit, the yield is independent of the initial conditions. This statement requires that the strong interaction states are 
in thermal equilibrium from the initial time U — > — oo. As discussed above, this does not apply to a QGP produced 
from a collision and thermalized at ~ 1 fm/c after the collision with a finite lifetime of a few fm/c. 

The detailed analysis of the dynamics of the virtual photon cloud highlights the difficulty and ambiguity in trying 
to separate the contribution from the virtual photon cloud from the contribution that grows in time during a finite 
lifetime. In the vacuum the virtual cloud relates the physical to the bare state, thus subtracting the vacuum virtual 
cloud amounts to studying the propagation of physical photons. However, the virtual cloud in the medium dresses 
a physical particle into a quasiparticle, hence the time dependent terms associated with the in-medium virtual cloud 
are actually describing the dynamics of formation of a plasmon quasiparticle. 

In section HVDI we analyzed the formation of the virtual photon cloud, and recognized that the negative frequency 
contribution features divergences associated with the virtual cloud of the vacuum as well as of the medium. Clearly 
these divergences are unobservable and must be subtracted, leaving solely finite contributions to the real time yield. 

The subtraction of the divergences associated with the virtual cloud is ambiguous during a finite time interval and 
we defined the real time yield by eq. (|IV.34f) to include only the positive and negative frequency contributions that 
are suppressed by the Bose-Einstein distribution function. Such definition leads to a finite and manifestly positive 
photon number density, but is not the only possible definition. 

In particular the subtraction i|IV.30|) of the virtual photon cloud of the vacuum, namely the T = contribution to 
the photon yield, entails that the quarks are asymptotic states in the infinite past and the photon cloud was built up 
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during the evolution of these asymptotic states from ti — > — oo up to the collision time. However, this assumption 
does not correspond to the actual QGP physics. 

Before the collision, quarks and gluons must be described as partons confined inside the nuclei and in terms of their 
distribution functions. Associated with the charged partons there is a photon cloud with a distribution function that 
depends on the charged parton distribution function in the nuclei. 

What happens to this photon cloud during and after the collision?, obviously this question is very difficult to address 
quantitatively, however the following are some possibilities 32]. 



i) The photon cloud is shaken-off by the collision resulting in a flash of photons during the pre-equilibrium stage. 
These photons have a different origin from those being emitted from parton-parton annihilation or bremsstrahhmg 
discussed in ref. 24]. If the photon cloud is shaken off after the collision, it will form again during the time between 
the collision and the thermalization of the plasma ~ lfm/c because of the electromagnetic interaction. 

ii) The partons de-confine during the collision becoming free and the photon cloud is not modified either by the 
deconfincment or by the parton-parton re-scattering that leads to a thermalized QGP. 

Clearly which of these (or other) possibilities describes the actual physics of the collision cannot be assessed with 
the current level of theoretical understanding. 

The importance of these questions cannot be underestimated. If the lifetime of the QGP were truly infinite, the 
photon yield at long times would be insensitive to the initial conditions and the contribution to the total yield from 
the initial stage will be of Q(l/t). However, given the short lifetime of the QGP, the initial condition not only is 
important but it bears an imprint in the spectrum. During the finite lifetime of the plasma there is no clear and 
unambiguous separation between the photons produced at a constant rate, those that are produced at a slower rate 
and those that are associated with the virtual photon cloud in the asymptotic long time limit. 

Although a detailed understanding of these issues is lacking, we can provide an approximate description that will 
model the essential ingredients. Thus, we now modify the choice of the initial density matrix in order to account for a 
period of electromagnetic dressing of the strong interaction eigenstates during the pre-equilibrium stage between the 
collision and the onset of thermalization. 

For this purpose, we now revisit the Gell-Mann-Low theorem[3l| that obtains the exact eigenstates of the total 
Hamiltonian from the free field in states in terms of the adiabatic M0ller wave operator. 

Define the M0ller wave operator 



where Hi(t) is the interaction Hamiltonian in the interaction picture of Hq. The Gell-Mann-Low theorem asserts that 
if the states |n) are eigenstates of the Hamiltonian Hq, then the states 



are eigenstates of the full Hamiltonian H. The M0ller wave operator adiabatically dresses the non-interacting bare 
state to be the full interacting dressed state during an infinite period of time. If quarks and gluons were truly 
asymptotic states the process of scattering, thermalization and photon emission would indeed be consistent with the 
dressing of the bare states from the infinite past. However, quarks and gluons are not asymptotic states, furthermore 
the actual collision involves bound states of quarks and gluons which are liberated after the nucleus- nucleus collision 
but the adiabatic hypothesis is not suitable to describe the process of formation or break up of bound states. Namely, 
the adiabatic hypothesis which is the basis of the S-matrix approach is not suitable to describe the dynamics of 
confinement and deconfinement. 

If partons are freed after the collision, the dressing process has to take place during the pre-equilibrium stage either 
completely, if in the process of the collision the partons shed their virtual cloud, or partially if the virtual cloud of 
photons that partons carried as bound states is also carried after the collision. This discussion brings to the fore 
the difficulty in separating unambiguously the virtual cloud from the observable photons during a finite lifetime, and 
manifestly makes clear the inadequacy of the S-matrix approach to describe any physical process in a QGP of a finite 
lifetime. 

Given that there is no current understanding of these issues, we now provide an approximate description of the 
dressing between the time of collision and the onset of local thermodynamic equilibrium. 

For this purpose, we define the interaction picture (with respect to Hq) quasi adiabatic interpolating states up to 
lowest order in the electromagnetic coupling as 
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The physical interpretation of these states is that the electromagnetic interaction dresses the eigenstates Hq during a 
time scale A natural time scale to describe the pre-equilibrum stage, between the collision and termalization is 
about 1 fm/c. 

These states interpolate between the exact in states when r — > + and the eigenstates of Ho for T — > +00. We 
note that these are entangled states in the sense that they are not simple tensor products of QCD and photon states, 
therefore these initial states are correlated. 

We can now construct an initial density matrix at time to m terms of these states that is reminiscent of a thermal 
density matrix but in terms of the dressed states (|V.3(1 . namely 



(*o) = 5>"^- K;7)K; 7 |- (V.4) 



The interpretation of this initial density matrix is that the QCD eigenstates have been dressed by the electromagnetic 
interaction on a time scale 1/r which describes the time between the collision and the onset of thermalization. 

It is important to note that this initial density matrix does not describe a state of thermal equilibrium under the 
strong interactions because it does not commute with Hqcd since the quark electromagnetic current that enters in 
the definition of the dressed states does not commute with Hq C d ■ 

We highlight this important point: any initial density matrix that includes the photon cloud, a result of the electro- 
magnetic interaction does not commute with Hqcd, hence it cannot describe a state in thermodynamic equilibrium 
under the strong interactions. The only manner to construct an initial density matrix with photons in the initial 
state and in thermal equilibrium under the strong interactions is for this density matrix to be factorized into a tensor 
product of a density matrix of pure QCD and a density matrix of free photons. The density matrix given by eq. 
(|iV.l|l is one such (the simplest) case. 

Once the initial density matrix IV. 41) is specified, its time evolution is completely determined by the full Hamiltonian 
and given by eq. (|III.17(I . The evolution of the number operator in time is therefore given by (|III.18jl . which upon 
inserting a complete set of eigenstates of Ho leads to the following result, 
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(V.5) 



The expression ea. l|V.5|) clearly coincides with ea. (|IV.10(l in the T — ► +00 limit. In the opposite limit r — * + it 
becomes 

dN r (t) u ^-oo 1 /ImII T (fc,^ = fc) ^ f)| r dcj ImU T (k,m) r 1 \ (y g) 



d 3 xd 3 k (2ir) 3 k \ e T - l J tt e* - 1 (<j - k) 

Thus, in the limit r — > + the yield agrees with ea. (|IV.12|) when the initial time to is taken to — 00 and the photo- 
production rate coincides with the result from the S-matrix calculation, as it must be. 
Furthermore, using the identity eq. I|IV.11|) . the long time limit t — to — * +00 yields, 

dN r (t) t-to^-oo 1 / lmn T (k,m = k) f f duo Imn T (fc,^) 1 | 

d^¥k Wfk{ TTTI L<_i ° + r l+ J - ef-1 V J^Wf ■ (VJ) 

Since is the time scale between the collision and the onset of thermalization, eq. (|V.7(I coincides with the S-matrix 
calculation in the long time limit from the instant of the collision. Therefore, this initial preparation is a physically 
acceptable description insofar as it reproduces the asymptotic long time limit. 

We note that the eq.JV3J) does not vanish at t = to because of the first term in the bracket, which is time 
independent. The value of the photon yield (|V.5Jl at t = to, namely the contribution determined by the first term in 
the bracket in eq. can be interpreted as the total number of photons (per unit phase space) created during the 

time scale These photons correspond to the virtual cloud as well as the observable photons emitted during the 
pre-equilibrium stage. This is precisely the physics that the density matrix in terms of the quasi adiabatic states is 
meant to describe. 

This interpretation becomes clear in the limit r — > in which case 

' T =° irT- 1 5{oj-k)+V . X M2 +0(r), (V.8) 
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the principal part leads to the divergences associated with the virtual cloud and since T _1 is the time scale of 
preparation, the term with the delta function gives the real photons produced during the time scale 

We then subtract the time independent term (first in the bracket) in eq. I)V.5(I . (which does not contribute to the 
rate). We thus define a photon number that vanishes at the initial time to and that is independent of the initial 
photon cloud, namely 
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The interpretation of this definition is gleaned from the expression 
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with ip—Ml given by eq. (|V.5|I . Thus, the subtracted number is obviously the photon yield between the time at which 
the plasma is thermalized t (after the collision) and the time t. Hence, this definition neglects the virtual photon 
cloud in the initial state and assumes that the observable photons produced during the pre-equilibrium stage leave 
the plasma. 

The extra powers of to — k in the denominator in eq. (|V.9fl render the total yield finite. 

The subtraction of the time independent term in eq. (|V.5|I has accounted for the divergent contributions of the 
virtual cloud of the vacuum and the medium, leaving a finite result for the time dependent yield. 

Therefore, this initial preparation and the subtraction of the total photon number at the time of thermalization, 
provide a possible systematic framework to approximate the physics of the initial state. 

In order to assess the (finite) contribution of the virtual cloud, it proves convenient again to separate the positive 
and negative frequency contributions. We then obtain 
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The term 



dN)?>{t) 



rfil ^ d3fc includes the vacuum contribution given by 7To(w, k) in ImIiT(fc,w). However, now this vacuum 
contribution is finite but of order V 2 /k 2 thus leading to a divergent number of photons which must be identified with 
a vacuum contribution to the virtual cloud and must therefore be subtracted. 
The remaining term 
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is finite and only depends on the medium. 

A lengthy but straightforward analysis of the time average (neglecting the oscillatory functions) yields the following 
result valid in the limit k 3> T,T 
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plus terms that are exponentially suppressed for k 3> T. Thus, the total yield, and the energy radiated in photons 
are finite. 
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FIG. 8: The contribution of eq. ifvTTKjl for k = 0.4 Gev/c (left panel) k = 3 Gev/c (right panel) for Y = 0.2 Gev ; T = 0.3 Gev. 



Fig. [S] shows the vacuum subtracted term ^S^far^ for two values of the momentum for a range of parameters 
expected at RHIC. 

Comparing fig. [S] with fig. [21 it becomes clear that the oscillations in fig. [S] are on much longer time scales, in 
particular it can be gleaned from fig. [5] that the shortest oscillation time scale is of 0(1/ k) and that there are longer 
time scales. While the results leading to fig. Elare sensitive to the cutoffs because the integrals diverge, leading to very 
rapid oscillations, the frequency integral in eq. I|V.15J| is finite an the integrand falls off fast. Hence, the oscillations 
are not sensitive to large frequencies and are on time scales which are of the order of the lifetime of the QGP. 

Thus, the contribution given by ea. (|V.15|l is: i) finite and leads to a finite photon number and energy, ii) the real 
time dynamics is on time scales of the order of the QGP, certainly at least for < k < 2 — 3 Gev/c. Therefore there is 
a priori no reason to subtract this term from the yield and it must be considered on equal footing as the contribution 
from eq. I|V.12|) . 

Thus, the final expression for the photon yield with initial preparation on a time scale T -1 and after subtracting 
the virtual cloud and the pre-equilibrium yield is given by 

dN F {t) _ 1 [°° du 



J — | ImII r (A;,u;) n(uj) T + [u,k,t ■ 



d 3 xd 3 k (2tt) 3 k 

+ [Imn T (fc, uj) [1 + n(u)] - ImIT T (fc, to; T = 0)] T~ [w, k, t - t ] j , (V.17) 

where T [u>, k,t — to] are given by eq. (|VT4|l . This expression is one of the main results of this study. 

While we obtained this expression based on the analysis of the lowest order contribution, we advocate eq. I|V.17(I 
to lowest order in a em and all orders in a s as an effective description of the photoproduction yield during a finite 
time interval. This expression includes the initial state preparation and has the following important properties: 

• The divergences associated with the virtual photon cloud both in the vacuum and in the medium, as well as the 
photons produced during the initial stage prior to thermalization are subtracted. 

• The total yield as well as the energy are finite. 

• The limit t — to —> +oo correctly reduces to the photon production rate obtained from the S-matrix calculation. 

• The limit r — > also leads to the S-matrix result for the rate. This is expected because T — > corresponds 
to an infinitely long preparation stage, which implies an infinitely long time interval, for which the S-matrix 
calculation applies. 
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• The initial preparation time scale T _1 is a parameter that describes in an effective manner the time scale between 
the collision and the onset of the QGP in LTE. It can be used as a fitting parameter for phenomenological 
purposes. 

Figure [5] displays the yield given by eq. (|V.17|) for several values of the initial preparation time scale for values 
of k, T for which the HTL approximation is valid. As mentioned before the HTL approximation gives the leading 
contribution for k <C T and does not depend on the vacuum contribution. 

The values T — 0.2 Gev corresponds to a time scale of about 1 fm/c which describes the time scale between the 
collision and the onset of a thermalized QGP. It is clear from the figure that change in the yield is rather minor for 
long-wavelength photons even in the case of an extremely long preparation time scale. Thus for k <SC T the real time 
yield from the lowest order 0(a em oP s ) is of the same order as the S-matrix yield during the lifetime of the QGP and 
is rather insensitive to the preparation time scale. 




FIG. 9: Comparison between the real time yield with initial preparation given by ea. l)V.17(l for k — 0.1 Gev; T = 0.5 Gev as a 
function of t — to (in fm/c) in the HTL approximation eq. 11 V.17I . 

Figure lTUI shows the comparison between the real time yield with the full one-loop photon polarization [eqs. (jIV.201 - 
IIV.23(l ] and the S-matrix yield as a function of time for a preparation time scale 1 fm/c for k = 0.4, 3 Gev; T = 0.3 Gev. 
Comparing the right panel of fig. ^|to the case r = oo displayed in fig. 0]we see that they are qualitatively similar, 
with the quantitative difference in the overall scale. However, it is clear from these figures that the yield from the real 
time calculation from processes that do not contribute to the S-matrix rate, is of the same order of or larger than the 
yield obtained from the S-matrix expression during the lifetime of the QGP. 

To emphasize this point further for larger values of the momentum, fig. 1111 displays the real time yield for k = 
3 Gev/c ; T = 0.3 Gev for a wide range of the time scale for initial preparation. 

It is clear from this figure that while there are a few quantitative changes with respect to the case V — oo, 
qualitatively the results are similar and all of the same order. 

Figure IT2"1 shows the spectrum of photons produced during the lifetime of the QGP expected at RHIC ~ 10 fm/c 
for a preparation time scale of 1 fm/c and T = 0.3 Gev. The left panel compares the real time and S-matrix yields 
vs. k. The right panel displays the logarithm of the yield vs. the logarithm of wavevector for the real time case only 
and clearly displays the power law fall of ~ k~ 5 for large momenta (k ^> T,T) as predicted by eq. I|V.16|) . These 
figures suggest a crossover from an exponential to a power law fall in the spectrum of the real time yield, the crossover 
ocurring at a value k c which depends on T. We find numerically that for r ~ 0.2 Gev k c ~ 2.7 — 3 Gev/c resulting 
in a marked flattening of the spectrum. We also find numerically that k c decreases upon increasing T. The power 
law dominance is a telltale of the contribution of the term in the real time yield eq. (|V.17|) that does not feature a 
Bose-Einstein distribution function which leads to an exponential suppression. As discussed above this power law 
leads to a finite number of photons and a radiated energy. In the infinite time limit the contribution that leads to 
this power law would be identified with the finite distribution of photons in the virtual cloud, but as analyzed and 
discussed in detail above, during the finite lifetime this term cannot be separated from the other contributions and 
enters in the yield on the same footing. 
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T=0.3 Gev 




t-t (fm/c) t-t„(lm/c) 

FIG. 10: Comparison between the real time yield N r t(t) with initial preparation given by eas. ljvTFfl with T = 0.2 Gev ; T = 
0.3 Gev and the S-matrix yield as a function of t — to (in fm/c) for k — 0.4 Gev (left panel) and k = 3 Gev (right panel). 
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FIG. 11: Comparison between the real time yield N r t(t) with initial preparation given by eas. IIV.17ll for several values of T, for 
k = 3 Gev/c; T = 0.3 Gev. 

Thus, the power law spectrum at large momentum is a hallmark of the processes that contribute during the finite 
lifetime of the QGP and that cannot be captured by the S-matrix approach. 

We have also studied the energy [see eqs. (UV.41|) -( |IV.42|l ]. which results in expressions similar to those given by eq. 
(|IV.45p but the positive and negative frequency contributions are replaced by those in the subtracted yield (|V.17|> . 
The numerical study of the energy reveals minor quantitative changes with respect to the results shown in figs. IBI7I 
The contributions of the terms which are not exponentially suppressed by the Bose Einstein distribution function 
begin to become important when the momentum is of order k ~ 3 — 4 Gev/c at which point all contributions are very 
small. The momentum integrals that lead to the energies are dominated by momenta < 1.5 — 2 Gev/c for T ~ 0.2 Gev. 
Consequently, the results for the energies from the initial density matrix with the initial stage preparation are very 
similar to the results displayed in figs. 16171 with an overall small change in the scale. 
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2 4 6 8 10 1 2 3 4 

k(Gev/c) Ln(k/(Gev/c)) 

FIG. 12: Comparison between the logarithms of the real time yield lniV r t(fc, t) with initial preparation given by eas. l|V.90 with 
F = 0.2 Gev and the S-matrix yield In NsM{k,t) as a function of k (in Gev/c) for T = 0.3 Gev. 

VI. CONCLUSIONS, DISCUSSIONS AND IMPLICATIONS 

In this article we studied the direct photon production from a QGP in local thermodynamic equilibrium during 
a finite lifetime. After discussing the shortcomings of the usual approach based on the S-matrix calculation of the 
emission rate, we focused on describing photon production directly from the real time evolution of an initial density 
matrix. 

The main premise of this study is that there are processes that contribute to the direct photon yield during the finite 
lifetime of the QGP but that are not captured by the S-matrix approach. We highlighted this point by restricting our 
study to the lowest order contribution of 0(a em ) to the yield. While this contribution is subleading in the asymptotic 
long time limit (that is if the QGP were truly a stationary state of infinite lifetime) it does contribute to the yield 
during a finite lifetime. 

We began our study by first considering an initial density matrix that describes a QGP in LTE with no photons 
in the initial state, compatible with all of the assumptions in the literature that lead to the S-matrix calculations of 
the emission rate. This study revealed the important aspect of the dynamics of formation and build-up of the virtual 
photon cloud. We highlighted that the finite lifetime of the transient QGP results in that the photon spectrum retains 
information of the initial state and pointed out the inherent ambiguities associated with the separation of the virtual 
and observable photons during the short timescale between formation and hadronization. 

While the terms that yield divergences in the number of photons and energy from the virtual photon cloud of 
the vacuum and in the medium can be identified, there is no unambiguous manner to subtract these from the finite 
contributions during a finite lifetime. Within this choice of initial state, we defined the yield subtracting the divergent 
contributions associated with the virtual clouds but also finite, time dependent terms because there is no unambiguous 
manner to extract these. The resulting yield clearly shows that there are contributions from processes that cannot be 
captured by the S-matrix approach, but that contribute to the direct photon yield during the finite lifetime. 

In particular, within the assumption of no initial photons, our study revealed that even after subtracting finite 
contributions, the yield from lowest order processes is of the same order of or larger than those obtained in the 
S-matrix calculations. 

We then provided an effective description of the preparation of the initial state by constructing the initial density 
matrix of a thermalized QGP in terms of quasi-adiabatic states obtained by electromagnetic dressing the QCD 
eigenstates over a time scale This time scale describes in an effective manner the pre-equilibrium stage between 
the deconfinement of partons and the onset of a thermalized QGP, it is expected to be of the order of 1 fm/c. This 
quasi adiabatic initial condition on the density matrix allows to extract the (divergent) contribution from the virtual 
cloud built-up as well as the observable photons emitted during the pre-equilibrium stage. The main result of this 
analysis is eq. (|V.17|) which provides an effective but systematic formulation of direct photon production from a 
QGP in local thermal equilibrium directly in real time. There are many advantages of this formulation over the usual 
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S-matrix approach: 

i) it docs not suffer from the caveats associated with the S-matrix approach described in section ITU. 

ii) it describes photon production in real time as an initial value problem, consistently with hydrodynamics [2fij . 

iii) it includes a description of the initial state in an effective manner in terms of the parameter T which is associated 
with the inverse of the time scale between the collision and that of thcrmalization. 

This is the time scale during which partons are almost free and parton-parton scattering brings the deconfincd 
partons to a state of local thermodynamic equilibrium. 

The final expression for the real time photon yield given by ea. l|V.17|) has several important properties: 

• It reproduces the result of the S-matrix result in the two limits in which it must be equivalent: the t — to — ► +00 
and the r — > limits. Both limits actually refer to a QGP of infinite lifetime, for which the S-matrix calculation 
applies. 

• It allows to separate the divergent contributions from the vacuum as well as the in-medium virtual cloud, along 
with the observable photons produced during the pre-equilibrium stage. Therefore this expression leads to a 
finite photon yield and a total finite radiated energy. 

• This description parametrizes the initial state in terms of a time scale which is a phenomenological pa- 
rameter. Namely it provides an effective description of the physics between the time at which nucleus-nucleus 
collision results in deconfined quarks and gluons and the time at which a QGP in LTE emerges. 

Phenomenological consequences: The main results of this study point out that during the finite lifetime of a 
QGP in LTE expected at RHIC or LHC there are processes that contribute to the yield that are not captured by the 
usual S-matrix approach. While we have focused on the lowest order such term, there are already many important 
consequences of phenomenological relevance: 

• The correct direct photon yield is actually larger than that calculated with the S-matrix approach. This is a 
consequence of the processes that contribute even at lowest order, and also processes missed by the S-matrix 
that arise from the region of u) 7^ k in the imaginary part of the photon polarization. Thus, the correct photon 
yield will be larger than the current estimates. A reliable estimate of the correction calls for a re-calculation of 
the imaginary part of the photon polarization for all u> 7^ k up to O (^a s In ^-^ . Such calculation is currently 
not available. 

• An important telltale of the processes that contribute to the yield during a finite lifetime is a power law spectrum 
of the yield of the form k . The coefficient of the power law bears information on the temperature as well 
as the time scale for thcrmalization L -1 of the plasma. This telltale is in striking contrast with the S-matrix 
yield which features an exponential fall-off. Depending on the value of the time scale T _1 this power law sets 
in for k c > 3 Gev/c with k c ~ 3 Gev/c for ~ 1 fm/c. While this power law spectrum may be an important 
signature, it sets in for a region of momenta in which a description of the QGP in LTE may break down. As 
mentioned above, the current data on the elliptic flow parameter i^fpT^fbH l23j| reveals large departures from 
hydrodynamics (+pQCD) which relies on a QGP in LTE for px > 2 Gev/c. In the region of momentum up 
to ~ 2 Gev/c the real time yield is almost indistinguishable from an exponential fall off, but begins to flatten 
towards the power law at about ~ 3 Gev/c. It is possible that the excess of photons and the flattening of the 
spectra in the WA98 data^j| may be explained by the processes studied here and that originate in the finite 
lifetime of a QGP. 

More questions: 

Our study indicates that direct photons from a QGP in LTE may not be a clean signature of the formation and 
evolution of the plasma as originally envisaged. The short transient nature of the QGP entails that the spectrum 
carries information on the initial, pre-equilibrium stage. However the electromagnetic properties of the initial state are 
largely unknown. The current estimates of the rate extracted from S-matrix calculations, which assume the existence 
of asymptotic states, and infinite QGP lifetime are not completely reliable, in particular, these are based on a weak 
coupling expansion in terms of a s but at the energy density conjectured to be achieved at RHIC a s ~ 0.24. Hence 
the estimates based on the S-matrix yield are at best qualitative. During the finite lifetime of the QGP processes that 
are completely neglected by the S-matrix approach give contributions to the yield that are of the same order as those 
of the equilibrium calculations, or even larger at large momenta, since the spectrum from the real time description 
features a power law fall off ~ fc -5 versus the exponential fall off of the equilibrium yield. 

Thus, in order to provide a phenomenologically reliable estimate the following questions would need to be addressed: 
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• What actually happens to the virtual cloud of photons after the collision??, are the virtual photons in the nuclei 
shaken off and if so does this result in a flash of photons during pre-equilibrium?. 

• The real time yield is sensitive to the structure of the photon polarization for u> ^ k 7 what is the full expression 
of the imaginary part of the photon polarization up to O ( a em a s In — ) ? 



• What is the range in momenta (fc-r) for which emission from a hydrodynamically expanding QGP is reliable?, 
if the data on elliptic flow for charged particles is extrapolated to photons (and in principle there is no reason 
to assume otherwise) then the local equilibrium description may only be valid up to fcx ~ 2 Gev/c. 

• We have discussed above that the terms that are identified with the virtual photon cloud in the medium 
asymptotically at long time are actually providing dynamical information on the formation of the quasiparticle 
in the medium. This is an aspect that has not been explored before, the formation of a quasiparticle in the 
medium does not require scattering and is to lowest order is independent of the mean-free path. This can be 
understood simply from the fact that the HTL approximation does lead to a plasmon quasiparticle but without 
collisional damping. Namely, the (transverse) plasmon is a consequence of Landau damping and not of any 
on-shell scattering process associated with a collisional width or a mean free path. As the produced photon 
traverses the medium it must necessarily carry with it a polarization of the medium that dresses the photon 
into a quasiparticle. What happens to this induced polarization once the plasma hadronizes??, is the virtual 
photon cloud of the QGP released in a flash during the hadronization transition?. If so a power fall off ~ fc~ 3 
in the spectrum is an unavoidable consequence of the formation and later dissipation of the virtual cloud in the 
medium. 

• We have studied the consequences of the finite lifetime to lowest order in the perturbative expansion corre- 
sponding to the one loop polarization. However for T/T c between 1 — 3, lattice data clearly shows that the 
quark-gluon plasma is not free or weakly interacting. Thus the next step in the program will consider self-energy 
and vertex corrections mediated by gluons, namely to higher order in a s . The strategy will be to compute the 
photon polarization including higher order corrections in a s and input its imaginary part in the final equation 
(|V.17|I . We expect to report on this study soon. 

• As we discussed above in section JnJ (see the discussion under 'caveats') the elliptic flow data suggests the that 
the hydrodynamic description is not valid for (transverse) momenta fcx > 2 Gev/c. Thus the assumption of LTE 
upon which all calculations of direct photon production from a QGP hinge, including the real time formulation 
studied in this article, will not be warranted for large momenta. However, and perhaps more importantly, 
at large transverse momenta it is expected that prompt photons produced during the pre-equilibrium stage 
during the pQCD parton-parton scattering will provide a large contribution to the total photon yield. Thus, as 
stated in ref.Q, the interpretation of the photon spectrum for large kr cannot be unambiguous. The current 
understanding of prompt photon production during pre-equilibrium is based on parton cascade calculations |24j 
which invoke a transport description and includes collisions via pQCD parton scattering cross sections. Such 
approach implicitly (and explicitly in the collision term) relies on an S-matrix description of the parton-parton 
collisions and is thus subject to a similar criticism described in section (|TTJ> above. Therefore a reliable estimate 
of the photon yield for large momenta requires understanding the dynamics beyond LTE and providing a reliable 
estimate of prompt photons from the pre-equilibrium stage. 

The experimental importance of electromagnetic probes of the QGP warrants a deeper study and assessment of 
these questions and in our view a re-evaluation of the current theoretical status on hard probes. 
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APPENDIX A: PHOTON POLARIZATION TENSOR 



The retarded photon polarization tensor is given by 

fLj,rei (■£ X t ) 



ie 2 {[Ji(x, t), Jj(f, t')])Q(t - t') . 



(A.l) 
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Introducing a complete set of simultaneous eigenstates of Hqcd and the total momentum operator P following the 
steps described in section HV Al above, we find 



(J i (x,t)J j (x',t')) = J d 3 p dco e -i?-(*-*)+M*-t') , (A.2) 

(Mx',t')MS,t)) = / d 3 p diu e -iH*-S')+Mt-t') a <ft u) , (A.3) 



where a^(p,u>) is given by ea. (|IV.5| l and 

<j<{p,u) = e-^(n q \Ji(6,0)\m g ) {m q \ J,-(0, 0)|n 9 ) 5 3 (p~Pn q +P,n q ) 8{u-E nq + E mq ) = o>faw) . (A.4) 

n q ,m q 

Introducing the Fourier representation of 0(i — t') we find the photon polarization tensor to be given by 



n^ ret (x-x\t-t>) = / — ^° e-^(--) e ^o(t-*') n ij>et (^p ) , (A.5) 



rf3 P rf P0 „-ip-.(£-x" 

(2tt) 3 2tt 

r+oo a > (v uj) [l — e' 3 "] 
n ij>et (p,p ) = (2tt) 3 / du Ijl n ■ - J . (A.6) 



Therefore the result, 



(2tt) a>(p,p ) = ^T^l • ( AJ ) 

We note that ImHij, ret (p,po) is an odd function of po with ImTIjj jT -et (p> Po > 0) > therefore a^(p,po) > for all 
values of po . 
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